APPROXIMATE ACOUSTIC CLOAKING IN 
INHOMOGENEOUS ISOTROPIC SPACE 

LIU HONGYU 

Abstract. In this paper, we consider the approximate acoustic 
cloaking in inhomogeneous isotropic background space. By em- 
ploying transformation media, together with the use of a sound-soft 
layer lining right outside the cloaked region, we show that one can 
achieve the near-invisibility by the 'blow-up-a-small-region' con- 
struction. This is based on novel scattering estimates correspond- 
ing to small sound-soft obstacles located in isotropic space. One of 
the major novelties of our scattering estimates is that one cannot 
make use of the scaling argument in the setting of current study 
due to the simultaneous presence of asymptotically small obsta- 
cle components and regularly sized obstacle components, and one 
has to decouple the nonlinear scattering interaction between the 
small obstacle components and, the regular obstacle components 
together with the background medium. 



1. Introduction 

A region is said to be cloaked if its contents together with the cloak 
are indistinguishable from the background space to certain exterior de- 
tections. Blueprints for making objects invisible to electromagnetic 
waves were proposed by Pendry et al. [22] and Leonhardt [16] in 2006. 
In the case of electrostatics, the same idea was discussed by Green- 
leaf et al. [I0] in 2003. The key ingredient for those constructions is 
that optical parameters have transformation properties and could be 
push-forwarded to form new material parameters. The obtained ma- 
terials/media are called transformation media, which we shall further 
examine in the current work for approximate acoustic cloaking in in- 
homogeneous isotropic space. 

The transformation media proposed in [TUl [22] are rather singular. 
This poses much challenge to both theoretical analysis and practical 
fabrication. In order to avoid the singular structures, several regular- 
ized approximate cloaking schemes are proposed in [7J [T31 [EH [T7J [23] • 
The basic idea is to introduce regularization into the singular trans- 
formation underlying the ideal cloaking, and instead of the perfect 
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invisibility, one would consider the 'near-invisibility' depending on the 
regularizes The works [7J and [23] are based on truncation, whereas 
in [131 [HI [17] , the 'blow-up-a-point' transformation in [10l[22] is regu- 
larized to be the 'blow-up-a-small-region' transformation. The perfor- 
mances of both regularization schemes have been assessed for cloaking 
of acoustic waves to give successful near-invisibility effects. Particu- 
larly, in [13J, the authors show that in order to 'nearly-cloak' an ar- 
bitrary content, it is necessary to employ an absorbing ('lossy') layer 
lining right outside the cloaked region. Since otherwise, there exist 
cloaking-busting inclusions which defy any attempts of cloaking. If 
one lets the lossy parameter go to infinity, heuristically we would have 
the sound-soft material lining, which is the one considered in [17] . All 
the aforementioned studies for approximate acoustic cloaking are con- 
ducted in the homogeneous background space, and there is no result 
available in literature for the more general case when the background 
space is allowed to be inhomogeneous. Nonetheless, it is worthy noting 
that the result in [13] for approximate cloaking of conductivity equation 
could be readily extended to the case with inhomogeneous background 
space by making use of the estimates in [B] for small extreme conduc- 
tivities; whereas the result in [13] could also be readily extended to 
the case with inhomogeneous background space by making use of the 
low-frequency estimates in [3] for the reduced wave equation. One of 
the key points for those extensions is that the corresponding studies 
could be reduced to the scattering estimates due to uniformly small 
objects, and then by scaling arguments, the studies could be further 
reduced to those having been considered in [B] and [3]. 

In this work, we shall consider the approximate cloaking for acous- 
tic waves in a very general and practical setting when the background 
space is allowed to be inhomogeneous but isotropic. We are mainly 
interested in the practical case that in the inhomogeneous space, there 
are both target objects one intends to cloak and non-target objects 
being uncloaked. By transformation-optics-approach, we construct the 
approximate cloaking devices by the 'blow-up-a-small-region' scheme. 
In order to overcome the cloaking-busts due to resonance, we imple- 
ment the sound-soft lining right outside of the cloaked region. In as- 
sessing the cloaking performance, the study is shown to be reduced 
to the scattering estimate due to extended objects with both asymp- 
totically small obstacle components and regularly sized obstacle com- 
ponents being presented simultaneously in an inhomogeneous space. 
One cannot make use of the scaling arguments as mentioned earlier 
and has to decouple the nonlinear scattering interactions between the 
small obstacle components, the regular obstacle components and the 
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background medium. By nonlinear scattering interaction we mean that 
the interaction between the scattered wave fields due to different com- 
ponents in a scattering system is a nonlinear process. Similar case 
with the background space being homogeneous has been investigated 
in [17], where boundary integral equations method is used to decou- 
ple the scattering due to the small obstacle components and regular 
obstacle components. For the current study, the presence of the inho- 
mogeneous medium makes the corresponding arguments rather tech- 
nical. We derive a novel system of integral equations underlying the 
scattering problem, which combines the volume potential operator of 
Lippman-Schwinger type and single- and double- boundary layer po- 
tential operators. An asymptotical coupling parameter is also a crucial 
incorporation. Since the scaling arguments does not apply here, an- 
other major difficulty one need to handle is that the domain of the 
underlying PDE is always in change in the asymptotic analysis. By ex- 
tensive use of the potential operators theory, we show that the scatter- 
ing contribution from small obstacle components is also asymptotically 
small with respect to their sizes, which justifies the near-cloak. 

In this paper, we focus entirely on transformation-optics-approach 
in constructing cloaking devices. We refer to [HI El Ell ESI EH] for 
state-of-the-art surveys on the rapidly growing literature and many 
striking applications of transformation optics. But we would also like 
to mention in passing the other promising cloaking schemes including 
the one based on anomalous localized resonance [20], and another one 
based on special (object-dependent) coatings [I]. 

The rest of the paper is organized as follows. In Section 2, we give a 
brief discussion on inverse acoustic scattering and invisibility cloaking. 
In Section 3, we collect the basics on transformation optics and apply 
them to the construction of approximate cloaking devices. Sections 4 
is devoted to the scattering estimates due to extended objects and the 
proof of the main result on approximate cloaking. 

2. Acoustic scattering and invisibility cloaking 

We first fix some notations that shall be used throughout the rest of 
the paper. For two domains D and Q in M n (n=2,3), D <<= Q means 
that D C D C Vt. xd denotes the characteristic function of D. Also, 
B r is reserved for a central ball of radius r in ~R n . For two relations 
TZi and 7^2,^1^ 72.2 will refer to 7Z\ < cR-i with c a generic constant 
which may change in different inequalities but must be fixed and finite 
in a given relation. ~ 7?.2 means that we have both < 7?.2 and 
7^2 ^ TZ-i- We shall also let c(zi,z" 2 , ■ ■ ■ denote a generic constant 
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that depends on the ingredients i\, i%, . . . , i\. The notations of function 
spaces is the usual one, e.g., L 2 (f2) denotes the space of square inte- 
grable functions on Q, and H S (Q) denotes the Sobolev space of order 
s on Q. 

Let cr = (cr* jf )™ J - =1 G Sym n be a symmetric-matrix-valued function on 
M. n and bounded in the sense that, for some constants < cq, C < oo, 

coft < ee(x)Z < Coft (2.1) 

for all and £ G IR n . Also, we let t] G L°°(IR n ) be a scalar function 

with 

7]{x) > r]o > ViGR". (2.2) 

In acoustic scattering, a -1 and r\, respectively, represent the mass den- 
sity tensor and the bulk modulus of a regular acoustic medium. Starting 
from now on, we denote by {IR n ; cr, 77} an acoustic medium as described 
above. We shall assume that the inhomogeneity of the acoustic medium 
is compactly supported, namely, a = I and rj = 1 in IR n \B with B a 
bounded domain in MJ 1 . In R n , the time-harmonic acoustic wave prop- 
agation is governed by the heterogeneous Helmholtz equation 

n 

di{a ij d jU ) + u 2 r]u = f, (2.3) 

where supp /nB represents a compactly supported source/ sink. Here, 
uj > represents the frequency of the wave scattering. In the following, 
we shall write {B; a, 77, /} to denote the scattering object including the 
inhomogeneous medium and the source/sink supported in Q. We seek 
solutions admitting the following asymptotic development as \x\ — > 
+00 

iuj\x\ ( 1 ^ 

< x ) = elX< + N (n-l)/2 [ A V, 0', w) + O(j^) j , (2.4) 

where 9,9' G S n_1 and £ = u8'. A(8,8',uj) is known as the scattering 
amplitude, which depends on the direction 9' and frequency u of the in- 
cident wave u l := e tx '^, observation direction 9, and obviously, also the 
underlying scattering object {Q; a, 77, /}. In inverse scattering theory, 
one intends to recover the target object {B; a, T], /} by knowledge of 
A(9,9',lo) (cf. [Hd5]). In the sequel, we use A({B; a,r], f}) to denote 
the scattering amplitude corresponding to {B; a, r\, /}. In this context, 
an invisibility cloaking device is introduced as follows. 

Definition 2.1. For a given regular background/reference medium {B; 
a ,r] }, let Q and D be bounded domains such that D <e Q d B. 
Q\D and D represent, respectively, the cloaking region and the cloaked 
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region. {fi\.D; a c , r] c } is said to be an invisibility cloaking device for the 
region D with respect to the background space {B; cr , r] } if 

A({B; a e , Ve , f e }) = A({B; a , rj }), 

where the extended object {B; a e , r) e , f e } = {D; cr a , r] a , f a }®{Q\D; a c , r] c } 
©{B\f2; cr , r/o} with {D; a a , r] a , f a } arbitrary but regular. Here, we use 
© to concatenate separate scattering objects. 

According to Definition ^. 1[ the cloaking medium a c , r] c } makes 

the target object {D; cr a , r] a , f a } indistinguishable from the background 
space {B; cr ,?7 }, and thus invisible to the scattering detections. 

As can be seen from our subsequent discussion in Section [31 one has 
to implement singular cloaking medium in order to achieve the ideal 
cloaking (see also the related discussion in [13] and [?]). This poses 
many challenges to both mathematical analysis and physical realiza- 
tion. In order to construct practical nonsingular cloaking devices, it is 
natural to incorporate regularization by considering approximate cloak- 
ing. We conclude this section by introducing the notion of approximate 
acoustic cloaking. 

Definition 2.2. Let {B; a , 1] }, Q and D be given as in Definition ^. II 
Let p > denote a regularizer and e(p) be a positive function such that 

e(p) -> as p -> + . 

{Q\D; o~P,r)P} is said to be an approximate invisibility cloaking for the 
region D if 

|| A({B;a e , Ve J e })-A({B; a , Vo }) || = e(j>) as p -+ 0+, (2.5) 

where the extended object {B; a e , r/ e , f e } is defined similarly to the one 
in Definition [2J] by replacing {Q\D; a c ,r] c } with {Q\D;a^,rj^}. 

According to Definition I2.2[ with the cloaking device {Q\D; o~p,i]P} 
we shall have the 'near-invisibility' cloaking effect within e(p) depend- 
ing on the regularizer p. 

3. Transformation optics and cloaking by sound-soft 

layer lining 

The transformation-optics-approach for the construction of cloaking 
devices critically relies on the following transformation invariance of 
the Helmh- oltz equation. Let x = F(x) : Q — > Cl be a bi-Lipschitz and 
orientation-preserving mapping. For an acoustic medium {Q; a, rj}, we 
let the push-forwarded medium be 



{tt; a, fj} = F*{Q; a, r)} := {Q; F*a, F*i]}, (3.1) 
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where 

a(x) = F*a(x) := -jMa(x)M T \ x=F -i^ ^ 

f)(x) = F*r](x) := r}(x)/J\ x=F -i {S:) 

and M = (dxi/dxj)? )j=1 , J = det(M). Assume that u G H l {tt) is a 
solution to the Helmholtz equation associated with {Q; <j,t]}, namely, 

V • (ctVm) + co 2 7]u = on f2, 

then the pull-back field u = (F _1 )*u := u o F^ 1 G H 1 ^) verifies 

V ■ (5-Vtt) + cu 2 fju = on (l, 

where we use V and V to distinguish the differentiations respectively in 
x- and ^-coordinates. We refer to [131 HZ] for a proof of this invariance. 

Next, we shall give the construction of the cloaking device. We 
first fix the setting for our study. Let {M"; o"o,^o} be a regular back- 
ground/reference space. Throughout, we shall assume that the back- 
ground space medium is isotropic, i.e., <To is a multiple of a scalar func- 
tion and the identity matrix. In the following, we treat <To as a scalar 
function which should be clear in the context. Let (Jq{x) G C 3 (lR n ) and 
?7o G C 1 (R n ), and supp (1 - cr ), supp (1 - r) ) <s B. Let £>, G <e B be 
C 2 -domains in IR n such that ~R n \(D U G) is connected and D is of the 
form D = U l k=1 D k , where each D k is simply connected. It is assumed 
that D k s are disjoint and separate from G. Let z k G D k be points such 
that 

\z k - z k i\ > d > 0, Vk^k', (3.3) 

and let and O k be C 2 -domains such that P k s and O^'s are simply 
connected and 

z k eP k <£ D k , G O fc . 

Define 

i 

U={JU k , U k = z k + O k . (3.4) 
fe=i 

Here, we assume that M. n \(U U G) is connected and 

distn«(E7, G) >r > 0, dist R n(C/ jfe , > f > 0, V fc ^ fc'. (3.5) 

For a sufficiently small p > 0, we let 

i 

u p = U ^ = 2fc + ' o0fc - ( 3 - 6 ) 

k=l 

Clearly, the parameter p determines the relative size of U kjP . Let p 
be sufficiently small such that U ktP d P k . We suppose that there 
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are orientation-preserving C 2 -diffeomorphisms F k : D k \U ktP —> Dk\Pk 
such that 

D k \P k = F k °(D k \U k , p ), F p k \ dDk = Identity, fc = l,2,...,Z. (3.7) 

One sees that blows up the small region U k , p to P k within D k . Here, 
we would like to mention that one of such mappings is given by 

Tix) — (a + 61x1)7-7, a = r 2 , b = , r 2 > r\ > 0, 

\x\ r 2 - p r 2 - p 

which blows up the central ball of radius p to the central ball of radius 
r\ within the central ball of radius r 2 . 

We are ready to construct the approximate cloaking device. Let 
P k , D k \P k denote respectively the cloaked and cloaking regions. The 
cloaking media are given by 

{D k \P k ; a k>c , r) k , c } = (F£)*{D k \U k>p ; a , Vo }, 1 < k < I. (3.8) 

That is, we blow up a small region U k;P within D k , and derive the 
cloaking medium in D k \P k by pushing-forward the background medium 
in D k \U k ^ p . Next, one can put the target objects in P k , namely, the 
passive regular media or radiating sources. The last ingredient for our 
cloaking construction is to introduce a thin sound-soft layer lining right 
outside P k . Then, P k can be regarded as a sound-soft obstacle. Here, 
by a sound-soft obstacle, we mean a scattering object which prevents 
acoustic wave penetrating inside (similarly, the wave propagation inside 
P k cannot penetrate outside and completely trapped inside), and the 
wave pressure vanishes on its boundary; that is, one would the following 
homogeneous Dirichlet boundary condition 

u — on dP k . 

For such construction, the target object together with its cloaking is 

1 

o~k,c, Vk,c} © P, (3.9) 

fc=i 

where P = U l k=1 P k is the union of all target sound-soft obstacle com- 
ponents; and the extended scattering object is given by 

£ = S®G®{B\{DUG);a , Vo } (3.10) 

The near-invisibility of our construction is given in the following 
theorem. 

Theorem 3.1. For any fixed £ e R n , we have as p — > + , 

\\A(6, £; S) - A(6, £; G © {B\G; a,M) lU^) = 0(e(p)), (3.11) 
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where 



<p) 




(lnp) 



-l 



n = 3 



n = 2 



(3.12) 



Theorem 13.11 indicates that the cloaked components P shall be nearly 
cloaked, whereas the uncloaked component G remains unaffected, even 
though there is scattering interactions among different components. 

We shall prove Theorem 13.11 in the next section, which is a conse- 
quence of the following theorem. 

Theorem 3.2. Let U p and G be sound-soft obstacles and 



\\A(6, £; T) - A{6, £; G © {B\G; a , Vo}) = 0(e(p)), (3.14) 



where e(p) is given in A3.12\) . 

Theorem 13.21 indicates that for scattering due to obstacles in the 
space {B; a ,r] }, the scattering contribution from small obstacle com- 
ponents is also small in terms of their sizes. The proof of Theorem 13.21 
will also be given in the next section. To our best knowledge, the re- 
sult in Theorem 13.21 are completely new in literature. We would like 
to remark that the scattering estimates due to small inclusions have 
received extensive mathematical interests in literature (see, e.g. [2]). 
So, our result in Theorem 13.21 together with the techniques developed 
for its proof is of mathematical importance and significance for its own 
sake. 

4. Scattering estimates of small obstacles in isotropic 



We first prove Theorem 13.21 The scattering problem corresponding 
to T in (13.131) is given by 

| div(a Vw) + u 2 7] u = in R n \(U p UG), 

I u\d Up yjOG = 0, 



T := U p © G © {B\(U P U G); a , r) }. 
For any fixed £ 6 8™, we have as p — > + 



(3.13) 



SPACE 



whereas the one without the small inclusions U p is 




(4.2) 
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Let u s (x) := u(x) — e tx '^ and u s (x) := u(x) — e tx '^ be the so-called 
scattered wave fields, and w s denote either u s or u s . Then, it is well 
known that w s satisfies the Sommerfeld radiation condition (cf. [_j) 

lim |x| (n ~ 1)/2 {^ -ujw'} = 0. (4.3) 
Ix^oo d\x\ 

We refer to [T_, __] for related study on the unique existence of u G 
H}jSL n \{U p UG)) to gU)-g3D, and u G H} oc (W\G) to (jOti- flO]) . 
However, for our subsequent discussion, we shall derive integral rep- 
resentations to u and u. To that end, we first recall the following 
celebrated gauge transformation (see, e.g. |24j). Let L a u = div(crVw). 
One has 

a' 1 / 2 o L a o a' 1 ' 2 = A - 7, 7 = a^Aa 1 / 2 . (4.4) 

1/2 ______ ________ 

Setting v = a u and using (I4.4p . ( 14. ip and ( 14. 3 p are transformed into 
(A + oo 2 q)v = in M. n \(U p U G), 

^|<9C/ P U<9G = 0, 



where 

-1 -2 -1/2 a 1/2 (A a \ 

q = cr iio-uj 70, 70 = ^0 Acr o • ( 4 - 6 ) 

Moreover, v s := v — e lx '^ satisfies the radiating condition (14. 3p . Obvi- 
ously, q G C 1 (R n ) and 1 — q is compactly supported in B. By intro- 

~ 1/2 ~ ~ £ 

ducing f = cr m and v s = v — e tx '^, we would have a similar system 
for v and v 8 , which should be clear in the context. 

We shall make essential use of integral representations to v and v 
for our arguments. To that end, we briefly introduce some potential 
operators and fix some notations for our subsequent discussion. Let 
fl be a bounded C 2 -domain in M n . We write Q~ as Q and Q + as its 
complementary, and II = dQ + = dQ~ . The one-sided trace operators 
for Q + and fl~ are denoted by 7 + and 7", respectively. The normal 
derivative of a function u G // 1 (fi ± ) is understood in the usual way as 
£>^ = YTj=i v jl ± {^3 v )i where the unit normal v points out of Q~ and 
into fl + . We shall drop the + or — superscript if the two-sided traces 
coincide. The adjoint operators 7* and £>* are respectively defined by 

(7**0,0) = (^,70)n for G ^(E n ) and ip G H 6 ' 1 ^), < e < 2, 

(4.7) 

and 

(£?*V>, 0) = OA, ^0)n for G <?(R n ) and i/j G #(11). (4.8) 
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Let 

e Mx-y\ j 

$(x,y) = —. r for n = 3; -H^\uj\x-y\) for n = 2, (4.9) 

with iip (t) the zeroth order Hankel function of the first kind be the 
fundamental solutions to —A — u 2 , respectively, in M 3 and 1R 2 . Now, 
we define 

Qu(x)= I <5>(x,y)u(y)dy for x G M n , (4.10) 

and the single-layer potential SL and the double-layer potential by 

SL = £ 7 * and DL — QB*, (4.11) 
whose integral representations are given by 

SLi/;(x) = / $(x,y)ip(y) dy, 



In the sequel, for notational convenience, we shall write for a set W C 

R n , 

SL U)W ip(x) := / $(x,y)if>(y) dy for y G W, 
Ju 

and similarly for DLu,w- Let 

Sji := SLu,u, Dji := DL u ,u 

which are understood in the sense of improper integrals. We shall also 
need the volume potentials 



JCiP(x) := / ${x,y)[l-q{y)]il>(y) dy, (4.13) 

'B\(U P UG) 



K<f>(x) := / $(x,y)[l - q(y)}<p(y) dy, (4.14) 
Jb\g 

and 

JC{x):= f *(x,y)[q{y)-l]C{y)dy (4.15) 
■J u p 

We refer to [H [19] for mapping properties of the potential operators 
introduced here. Finally, we let 

n = B\(U p UG), £l = B\G, 
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and 



<9G:=E; r fc . 



p •' 



dU> 



k,p, 
I 



P 



k=l 



dU p - 



T k := dU k , r := |J T k = dU. 

k=i 

Lemma 4.1. Let v l (x) := e lx '^. The solution v G C 2 
i n \(U p U G)) to ggp awen 6y 



=f*(x) — w 2 

<9z/(n) 
<9$(x, y) 



+ 



\(U P UG) 

-i$(x,y) 



n \(u p u G)) n 

[i - ?(2/)My) dy 
V»i(y) ^ (4.16) 



— y) 



r P L du (y) 

where x G R n \(f/ p UG) and k = [e(p)]- 1 . /fere, G C{Q), ^ G C(E) 
and -02 ^ C(rp) are uniquely determined by the following system of 
integral equations, 

v + u 2 JCv - (DLz q - iSL^)^ x 

- ( DL r p ,n - inSLr si)^2 =P into, 



^2 



+ (DL Vpt v - mSX rpiS )V>2 = <?i on E, 



(4.17) 
(4.18) 

(4.19) 



+ (£>r p - iKS Tp )i/j2 = q2 on T 



p- 



where p(x) = v l (x) for x G 0, gi(x) = — p(x) /or x G E and, 02 (x) = 
— /or x G T p . Tae system ^.11^ - ^.19\) is uniquely solvable in 
C{Q) x C(E) x C(r p ) /or every p G (7(0), pi G C(E) and p 2 G C(r p ). 
It is also uniquely solvable in L 2 (Q) xC(S) xC(T p ) for every p G L 2 (Q) 
and, pi G (7(E) and p 2 G (7(r p ). 

We remark that a similar scattering problem to (14. 5|) is considered 
in [18] , where the combination of volume potential and boundary layer 
potentials are implemented to represent the scattered wave field. How- 
ever, the main concern of our study is the scattering behaviors due 
to the asymptotically small obstacles U p . The novelty of the integral 
representation (14.161) is the introduction of the asymptotically coupling 
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parameter k for layer potentials on the boundaries of the small obsta- 
cle components, which shall be crucial for our subsequent scattering 
estimates. The proof of Lemma 14.11 follows from a similar manner to 
that in [18] , which for completeness we briefly include in the following. 



Proof of Lemma \4-l\ Using the mapping properties of potential oper- 



ators (cf. [1]), along with the fact that (— A 4 



)$(z,y) 



S x , it is 



straightforward to show that for v given in f )4.16p . v G C 2 (M. n \(U p U 
G)) n C(R n \(U p U G)) satisfies (A + u 2 q)v = 0. Whereas by the jump 
properties (cf. [1]), (14. 18}) implies v |s = and (14.191) implies v\r p = 0. 
The radiation condition for v s is a direct consequence of the integral 
kernel <$>(x,y). So, we only need to show the well-posedness of the 
system (I4.17p - (l4.19p . which could be written as 



(A + K)u = p 



(4.20) 



where 





" I 




-DL Fp , n ' 




( v \ 




( V 


A = 





I 





, u= | 


1>1 


1—1 


2gi 










I 








V 2g 2 



(4.21) 



and 



K 



o; 2 /C iSL-^Q 
-2w 2 7 /C 2{Dj: - iSx) 
-2u; 2 7/C 2(DL^ Tp _ lSL v 



It is verified directly that 



I DLx >n 
/ 




2{DL Tp ^ - iKSLr p ,x) 
2(D Fp - ikS Tp ) 



DL T 

I 



(4.22) 



and K is compact in both L 2 (tt) x C(E) x C(T p ) and C(Q) x C(E) x 
C(r p ). So, A + K is an Fredholm operator of index 0. We only need to 
prove the uniqueness of the system (14.201) . and it suffices to show this 
in L 2 {VL) x C(S) x C(Y p ). Set p = 0. By the uniqueness of solution to 
the scattering system (14.51) . we see v = in Q. Let 



w ix) 



d&(x, y) 
dv{y) 
d$(x, y) 
8v(y) 



m$(x, y) 



V>i(y) d y 

Mv)dy, xer\(sug. 
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Then, by 1 14. 17ft . we see w|n = v\n = 0. By the jump properties of layer 
potential operators, we have 

_ dw 

- 7 w = ^i, -7 (— ) = i^i on S; 
_ <9w 

- 7 w = V> 2 , -7 (-^j) = i«^2 on T p ; 
By Green's formula, 



i ■> _ _ <9w 

i / IV'il ds = / (7 w )7 (7^7) 



|Vw| 2 — w 2 |w| 2 (ix 



G 



which implies ^1 = 0. Similarly, one can show that ^2 = 0. The proof 
is completed. 

□ 

Remark 4.2. Similar to Lemma l4TTj we know the solution v G C 2 (M. n \G)r\ 
C(W n \G) to the scattering problem without the small obstacle U p is 
given by 



V[X) =V (X — 04 



®(x,y)[i - q{y)]v{y) dy 



+ 



d$(x, y) 
du(y) 



-i$(x,y) 



(4.23) 



i)(y) dy xe R n \G, 



with uU G C(f2) and ^ G C(S) uniquely determined by the following 
system of integral equations, 



v + 04 Kv — (DL^ = — iSL^ = )^ = p in O 
- u 2 'y)Cv + (D s - iSs)^ = g on S, 



(4.24) 
(4.25) 



where p(x) = ?/(x) for x G fi, q(x) = —p(x) for x G E. For the 
subsequent use, we let 



V 




p 




, x = 





(4.26) 



Clearly, we have 

v = L _1 x, 

both in L 2 {Q) x C(E) and C(5) x C(£). 



(4.27) 
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We next derive the key lemma in proving Theorem I3.2L In the fol- 
lowing, the analysis is based on the space dimension being 3. Later, 
we shall indicate the necessary modifications for the two dimensional 
case. Henceforth, for an operator A : X i— > Y with X, Y being Banach 
spaces, we shall denote by ||A||£(x,y) the corresponding operator norm. 

Lemma 4.3. For p — > + , we have 

\\V - v\\ L 2 {n) = \\V - v\\ L 2 ( Ci\U p ) = °(p) ( 4 - 28 ) 

and 

\\1>i-1>\\cw = 0{p) and |Mc(i>) = 0(1). (4-29) 
Proof. The proof shall be proceeded in 6 steps. 

Step I. We introduce w = WiXsi\u + w iXu p G L 2 (Q) such that 

Wl =v in tt = Q\U P , (4.30) 

and 

w 2 (x)+oj 2 &(x,y)[l-q(y)]wi(y)dy 



JQ,\U P 

-{DL^ Up -iSL^ Up )^{x) (4 ' 31) 
-(DL rpiUp - iKSL Tp)Up )ip 2 (x) = p(x) x G U p . 

By (jOIj) . fTOOj) and ffl~TTI) . one can see that w G L 2 (Q) and w 2 G 
L 2 {U P ) satisfy the following system of operator equations 

w + u 2 JCw + uJ 2 ' Jw 2 —{pL^ q — iSLj, 

-( DL r p ,Q- iKSL r p ,n)' l l J -2=P infi, (4.32) 
w 2 + u 2 Jw 2 + u 2 K,w-(DL T , :Up - iSL^u^tpx 

-(DL Tp ,u p - iKSL rptUp )ip2 = V in U p . (4.33) 

Step II. Let Q(x) = x/p : U p U T p — > U U T. In the sequel, for x G T, 
we define 

(sMx) = Jm*,v)<Kv) dy, {d&)(x) = l^^y^y) d y> 

where §o(x, y) = 1/ (Att\x — y\) is the fundamental solution to —A. For 
G C(T p ), by using change of variables in integration, one has 

e w\x-y\ r e iwp\x'-y'\ 

-<t>{y) dy = p — t jr(j)(py') dy', 



x — y\ J r \x' — y' 

d ( t r /du(y)<j>(y) dy = d( 6 — ) / dv{y')<t>{py') dy', 

x ~ y\ / Jr \ \ x ~ y \ 



APPROXIMATE INVISIBILITY CLOAKING 



15 



where x' = x/p,y' = y/p G T. Using these along with power series 
expansion of exp{iup\x'—y'\}, one can verify directly that for G C(T p ) 

IK0- 1 )* o S rp o (e- 1 )^ - S° r o (e-^VIIocr) < pU\\c(r p ), (4-34) 

where in the above inequality we have identified <fi{x) G C(T P ) as = 
<P(px>) = (e- 1 )^ G C(r), and the fact that |M|c(r„) = || (e- 1 )*^^). 
Similarly, we have 

IKe- 1 )* o d Tp o (e- 1 )^ - d° t o (e-^Woir) < p 2 U\\c { t p) - (4.35) 

Next, by fl4?L9l . together with the use of fQ4j) and we have 



(e- 1 )*^ =[-/ + d° - is° + o^j-^e- 1 )* 



q 2 + u?^Kv 



(4.36) 



Here, we made use of the invertibility of (^I+D® — iS°) : C(T) i-> C(r) 
(cf. 0). 



Step III. Using f)4.36p . we consider 
{DL Tp ^ - inSL rp ^)ip 2 

=(£>Lr„ E - mSX r/ ,, s )e*{[~I + D° r - iS° + O(p)]- 1 



(4.37) 



(6 



We first note that 

dist]Rn(r, E) < dist R n(r p , £) < diam K n(B), 
and hence by direct verification 

||^D^r p ,s - iKSL rp p\\c(c(r p ),c(p)) = O(p) (4.38) 

and 

\\DLx iTp -iSLz,rMcw,CQ? p )) = 0(1)). (4.39) 
By letting v = vxn + ®Xu p £ £ 2 (fi), and using the fact that /C maps 
L 2 (Q) continuously into H 2 (Q) (cf. [1]), we also see that 



IMIc(n) < ll^llc(Q) ~ W^Wh^q) 



~ IHIz, 2 (n) 



\v\\v*(n) < \\w\\L2(ny 
By (O7HP011 . we have 

(DLr x - mSL Tp , s )^2 = A^v) + A 2 {^i) + A 3 {q 2 ), 



(4.40) 



(4.41) 
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where 



Ai(v) ={DL Tp ,x - tKSL rp ^e* ([±1 + D° r - iS° r 

x ~ (4.42) 



+ 0{p)]-\Q- 1 )\u 2 1 lCv) 
satisfying 

MiMllc(E)<p|kl|L 2 (n)<p|klL 2( ^; (4.43) 

and 

A 2 M =(DL Tp ^ - wSL Vp ^W + D°- iS° r + O(p)]- 1 

V v (4-44) 
(e- 1 )*(-DL Ei r p +z5 , L E , rp )ViJ 

satisfying 

\\M1>i)\\ <Pll^i||c( S) ; (4.45) 

and 

A 3 (q 2 )=(DLr p ,v-iKSL Vp>s )G* 

1 (4.46) 

(i-/+z)?-^ + o(p)r 1 (e- i rg 2 ) 

satisfying 

\\Mto)\\c<E) < p\\<h\\c(r p y (4-47) 
Hence, we see from (14.471) that 

WMMcp) < P- (4-48) 
Plugging (147361) into (147181) . and using (147371) - fOHD . we have 

[h+Dz-iS s + 0(p)]il> 1 -u 2 'y]Cv + A 1 (v) = q + 0{p) on S. (4.49) 

Step IV. In the sequel, we shall denote by DL° and SL° the double- 
and single-layer potentials with the integral kernel <&(x,y) replaced by 
&o(x,y). Also, we let 

Jo({x) = / <S> {x,y)[l-q{y)]({y) dy, 
Ju 

where q = (O -1 )*? G L°°(U). We know that / + Jo is bounded invert- 
ible from L 2 (U) to itself and (cf. [I]) 

IKJ + ^Jo)- 1 !!^^)^^)) < c(Z7,||g|| L =c (t0 ,w). (4.50) 
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Clearly, the bound in (I4.50p is independent of p. Now, using change of 
variables in integration and by a similar argument to that in Step II, 
we have from ( 14 . 3 3 1) 

(e- 1 )*^ = [/ + co 2 j + o{p)]-\Q- l y - u 2 Zw 

+(DLx tUp - iSLx,u p )ipi + {DL VpyUp - inSL rpiUp )tp 2 + V 

=Mi{w) + M 2 (ipi) + M 3 (ip 2 ) + Mi{p). 

We next assess the L 2 (C/)-norms of Ai\(w), A4 2 (tpi), ^3(^2) and 
Ai^p), respectively In (14.511) . 

Mi(w) = [I + io 2 J + 0( P )}- l (Q- l )*[-uj 2 1Cwl 
which satisfies 

IIMlHIIl-CeO < \\Kw\\ om < \\w\\ L ^ y (4.52) 

In (PT]) . 

= [I + oo 2 J + 0{p)]- 1 {Q- l y\{DL^ Up - z5L E)0p )Vi], 

Noting 

dist]Rn(E, r) < dist K n(S, U p ) < diam K n(B), 
one can easily show that 

\\(DL^ Up -iSL^u^Wau,) < ll^illc(s), 

and hence 

\\M 2 M\\l*(U) £ HV'illc(B). (4.53) 

For M.?,(ip 2 ) ) we have 

M 3 M = [I + uj 2 J + 0{p)}- l {®- 1 y\{DL Vp , Up - iKSLr ptUp )fo). 
It is verified directly that 

f / iui\x-y\\ 

{DL Vp , Up ^){x) = J^d yy^J /dv{y)Mv) d y 

/ iwp\x'-y'\ \ 

d(— f -)/du{y')Upy')dy\ 

r \W-V\J 
r e M x ~v\ r e ^pW-y'\ 

(SLr uM(x) = / 1 My) dy = p — jr^ipy') dy', 

Jr p \x-y\ Jr F V I 

where x G T p , y G U p and x' G T, y' G U . Hence, we have 

[DL Tp , Up - ikSL Tp ^ 2 = [DL° TtU + iSL° TtU + 0{p)]{Q- l y^ 
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and from which we further have 

\\M 3 (i> 2 )\\L*(u) < \\H\c(r p y (4.54) 

Finally, by 

Mi® = [i + tu 2 j + o{ P )}- l {Q- l y P , 

we see 

\\M&)\\ L * {U) = 0(1). (4.55) 
Now, by (l4T5il - fj435l) . we have 



(4.57) 



\\w2\\l*{u p ) = P 3/2 ||(© 1 )*^ 2 ||l2 (C7) < p* ,2 \\\Mi{w)\\ L 2 {u) 

+ \\M 2 (?l>i)\\vV) + WMsMWlhu) + \\M 4 (p)\\ L 2 {u ^ (4 ' 56) 
< P 3/2 \H\mn) + P 3/2 |l^illc (S ) +p 3/2 ||^||c ( r p ) + 0(p 3/2 )- 

Step V. We now consider the equation ( 14.321) . Let w = W2Xu p + ®Xn 
L 2 {Q). Using (14361) . we have 

\\<? w 2\\m(ci) = II -^llfla(fi) ~ IHI^^ = \M\lz(u p ) 
< P 3/2 \H\l^) +P 3/2 |l^illc;(E) +p 3/2 ||^||c ( r p ) + 0(p 3/2 ). 
By gSED, (ESD and (OOl) . one can show 
IWIcf,) = ||(0 _1 )>2||c(r) 
^ ll?2||c(r p ) + 11^11x3(0) + ll^illc(E) (4.58) 
^ ll?llc(n) + IMIl^) + ll^illc(s) + O(p), 
which together with (14.561) and (14.57ft implies 

\\w*\\i?<p p ) < P 3/2 |kllL 2 (n)+P 3/2 |l^i||c(E) + 0(p 3/2 ), (4.59) 

and 

ll^llfl>(fi) S P 3/2 lkllL 2 (n)+P 3/2 |l^illc (S ) + 0(p 3/2 ). (4.60) 

Let $ G ^(M 3 ) be a cut-off function satisfying •& = 1 on B and = 
on D with D a neighborhood of B. For <fi G 3t(M. n ), we have 

(DL r(3iR n'02,#)L2(Rn) = (^S*'02,#)L2(Rn) 

= (^ 2 ,^(^^)0) L2(rp) < ||^||^ ( rJ£,(^MI^(r,) 



< Pll^||c(r p )||(^)0lk 2 (D) 

< P||^2||c(r p )||^||L2(R"), 



(4.61) 



and from which one easily see that 

\\DL Tp ^2\\mn)<p\\H\c(r p ). (4.62) 
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Here, we have made use of the fact that Q* maps L 2 (D) continuously 
into H 2 (D). Similarly, we have 

(S£r p ,R^2,#)£ 2 (K«) = {G&y*i>2, #)i 2 (R") 

= (^ 2 , 7 (W)0)L 2 (r p ) < ||^2||L 2 (r p )||7(^)0||L 2( r p ) 
<P 2 ||^||c(r p) ||(^)0||c(r p ) (4.63) 
< P 2 ||^2||c ( r p )||(^)0||^(D) 
< p 2 ||^2||c(r p )||^0||L2(R"), 
and from which one has 



no 



(4.66) 



(4.67) 



\\ K S L r p ,M\L*(n) =p 1 \\SL rp ^ 2 \\ I p i f (i < p\\ip2\\c(r p )- (4.64) 
By (gS2) and (Oil) , we see that 

IIP L r p ,n - ^^.fl)^ 11^(6) < p||^||c(T P ), ( 4 -65) 
which together with (I4.58P further gives 

\\( DL r p ,n ~ m5L r p ,n)^2||L2 ( ^ } 
£ pIMIl2(q) +Pll^i||c;(s) + O(p). 
Obviously, (14.601) implies that 

ll^all < \\Jw 2 \\ 

<P 3/2 lkll^) + p 3/2 ll^llc (S ) + 0(p 3/2 ). 

Finally, (14321) . (141)61) and (ES3), together with the fact that q = q in 
fi, yield that {w,ipi) G L 2 (fi) x (7(E) satisfies 

[/ + c 2 /C + 0(p)]w - [DLzfi - iSL^ + 0{p)]^ 

(4.68) 

=p + O(p) in Q. 

Step VI. In the following, we further assess (I4.49p . We first note 

■yJCv = ■yJCw + ^Jw2 on E. (4.69) 
By (14.571) and (14.581) . one can show 

\\lJw 2 \\ C (Z) = || -T^llc(E) < \\fcw\\H*(a) , , 

<p 3/2 |kll^)+P 3/2 |l^lllc(E)+0(p 3/2 ), 

which together with (I4.69p . (14.491) and (I4.43P implies 
[^I+D^-iS^+0(p)]ip 1 -[u 2 jK:+0(j>)]w = q+0(j^ on E. (4.71) 
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Now, by (Qgjl and fl4\7H we see that v := G L 2 (fi) x C(E) 

satisfies 

[L + 0(p)]v = x + 0(p), (4.72) 
which by Remark 14.21 gives 

v = v + C(p). (4.73) 

That is, 

\\w-v\\ LH ^ = 0(p) and W l -$\\cp) = 0(p), (4.74) 
which together with (I4.58P implies 

||^l|crcr p) = 0(l). (4.75) 
It is trivially pointed out that 

\\v - v\\v>(a) < \\w - v\\ L 2 { n) = 0{p). (4.76) 
The proof is completed. 

□ 



Remark 4.4. In two dimensions, we shall have as p — > + 

\\v-v\\ L 2 m = 0((\np)- 1 ) 

and 

Ui-ncV) = 0((\np)- 1 ) and \\^ 2 \\ c{Tp) = <D{1). 

The proof is very similar to that for the three dimensional case in 
Lemma I4.3[ and the only modification we shall need is that at certain 
point, we shall make use of the following asymptotic expansion for 
$(x,y), 



(4.77) 
(4.78) 



$(x,y) =-H^\u\x-y\) = —In ' 



i 

+ - 



2n \x — y\ 4 

In VO \ \x-yrhx -. ; 

2tt 2 2tt V \x-y\ 

as | a; — y\ — > 0. Here, E denotes the Euler's constant. 

Proof of Theorem \3.2[ By noting u s = v s in IR n \B and (14.161) . one has 
by direct asymptotic expansion that 



A(9,^T) = -uj 2 t I e-^y[l-q(y)]v(y) dy 
n 



+ T 

+r 



. dv(y) 
Q e - iuj Q-y 

du(y) 



ipi(y) d y 



(4.79) 



ine 



-iuiO-y 



ifo(y) dy, 
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where To = e 1 ^ /y/8iru} for n = 2; and To = j- for n — 3. Similarly, by 
■£t s = a Q l ^ 2 v s and (14.231) . one also has 

A(9,£;G®{B\G;a ,r) }) 

f e-*"v[l - q{y)]v{y) dy 
In 



U 2 T 



(4.80) 



+ r 



-iujQ-y 



dv{y) 

Subtracting (14.801) from (14.791) . we have 

£; T) - A(6, £; G © {B\G; <r , 770}) 



- a; r 
+ r 
+ r 
+ w 2 r 



-***[!- q(y)][v(y) - v(y)] dy 



de~ iuje ' y 
du(y) 

Q e -iw8-y 



ie 



-iw8-y 



-iuid-y 



foiv) dy 



(4.81) 



du(y) 

e — °-v[l-q(y)]v(y) dy 



Applying Lemma I4T31 and Remark l4.4l to estimating the terms in (I4.8ip . 
along with straightforward calculations, we have (13.141) . The proof is 
completed. □ 

We are in a position to prove Theorem 13.11 

Proof of Theorem VJ . 1\ Let u G Hj- oc (M. n \(PUG)) be the scattering wave 
field corresponding to S = ^[ =1 {D k \P k ; a kiC , rj ktC } ©P©Gffi{B\(Du 
G); o-Q,r] }, and u G Hl oc (M. n \G) be the scattering wave field corre- 
sponding to G © {B\G; CT , ?7o}. Define F : W n \(U p UG) -> R n \(P U G) 
by 

F fc;P on D k \U p , 
Identity otherwise. 



F :-- 



(4.82) 



Obviously, F is bi-Lipschitz and orientation-preserving. By the trans- 
formation invariance of the Helmhlotz equation as we discussed in Sec- 
tion [3j one sees that F*u G Hl oc (M. n \(U p U G)) is the scattering wave 
field corresponding to T = U p © G © {B\(U P U G); do) r] }. Since F is 
the identity outside D = U l k=1 D k , we know 



A(e,{;£) = A(e,£;T). 



(4.83) 
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According to our discussion in Lemma 14.11 and Remark 14.21 we know 
u e C 2 {R n \G) and F*u e C 2 {R n \{U p U G)) n C(R n \(U p U G)). Then, 
by Theorem 13.21 

\\A(0, £; T) - A(9, £; G © {B\G; <r , 770}) H^s™-*) 

=C(e(p)) as p -> 0+, 

which together with (I4.83P yields 

|| £; £) - £; G © {B; a , r^o}) U^-i, = 0{e{p)). 

The proof is completed. □ 
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